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Abstract: Recent perturbative studies show that in 4d non-commutative spaces, the 
Q^' trivial (classically stable) vacuum of gauge theories becomes unstable at the quantum 

level, unless one introduces sufficiently many fermionic degrees of freedom. This is due to 
a negative IR-singular term in the one-loop effective potential, which appears as a result 
of the UV/IR mixing. We study such a system non-perturbatively in the case of pure 
U(l) gauge theory in four dimensions, where two directions are non-commutative. Monte 
5^ I Carlo simulations are performed after mapping the regularized theory onto a U(A^) lattice 

gauge theory in d = 2. At intermediate coupling strength, we find a phase in which open 
Wilson lines acquire non-zero vacuum expectation values, which implies the spontaneous 
breakdown of translational invariance. In this phase, various physical quantities obey clear 
scaling behaviors in the continuum limit with a fixed non-commutativity parameter 6, 
which provides evidence for a possible continuum theory. The extent of the dynamically 
generated space in the non-commutative directions becomes finite in the above limit, and its 
dependence on 9 is evaluated explicitly. We also study the dispersion relation. In the weak 
coupling symmetric phase, it involves a negative IR-singular term, which is responsible 
for the observed phase transition. In the broken phase, it reveals the existence of the 
Nambu-Goldstone mode associated with the spontaneous symmetry breaking. 
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1. Introduction 

Non-commutative (NC) geometry [1,2] has been studied extensively as a modification of 
our notion of space-time at short distances, possibly due to effects of quantum gravity [3]. 
It has recently attracted much attention since gauge theories on a NC geometry have 
been shown to appear as a low energy limit of string theories with a background tensor 
field [4] . At the classical level, introducing non-commutativity to the space-time coordinates 
modifies the ultraviolet dynamics of field theories (the interaction becomes non-local at the 
scale of non-commutativity), but not the infrared properties. The latter changes at the 
quantum level, however, due to the so-called UV/IR mixing effect [5]. As a consequence, one 
cannot retrieve the corresponding commutative theory by sending the non-commutativity 
parameter to zero in general. 
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The UV/IR mixing poses a severe problem in the renormalization procedure within 
perturbation theory since a new type of IR divergences, in addition to the usual UV di- 
vergences [6], appears in non-planar diagrams [5]. The finite lattice formulation [7] based 
on twisted reduced models [8] and their new interpretation [9] in the context of NC ge- 
ometry, regularizes such divergences — not only the ordinary UV divergences but also the 
novel IR divergences. It therefore provides a non-perturbative framework to establish the 
existence of a consistent field theory on a NC geometry. In ref. [10] a simple field theory, 
2d U(l) gauge theory on a NC plane^, has been studied by Monte Carlo simulation, and 
the existence of a finite continuum limit has been confirmed. The ultraviolet dynamics is 
described by the commutative 2d U(oo) gauge theory. On the other hand, the infrared dy- 
namics would be described by the commutative 2d U(l) theory, were it not for the UV/IR 
mixing effects. The "dynamical Aharonov-Bohm effect" observed in this regime clearly 
demonstrates the non-equivalence to the commutative 2d U(l) theory. 

As an interesting physical consequence of the UV/IR mixing in the case of scalar field 
theory, ref. [13] predicted the existence of a "striped phase", in which non-zero Fourier 
modes of the scalar field acquire a vacuum expectation value and break the translational 
invariance spontaneously. This prediction, which was based on a self-consistent Hartree- 
Fock approximation, was supported later by another study using an effective action [14], 
and discussed also in the framework of a renormalization group analysis in 4 — e dimen- 
sions [15]. Finally, the existence of this new phase was fully established by Monte Carlo 
simulations [16-20]. In ref. [19], in particular, the explicit results for the phase diagram and 
the dispersion relation were presented in d = 3. The dispersion relation in the disordered 
phase shows a positive IR singularity due to the UV/IR mixing, which shifts the energy 
minimum to a non- vanishing momentum. As one changes the mass parameter in the action, 
the mode at the energy minimum condenses, which yields the corresponding stripe pattern. 
The existence of a sensible continuum limit is suggested by the scaling of various physical 
quantities. In particular the average width of the stripes stabilizes at a finite value in the 
continuum limit. In even dimensions with the non-commutativity tensor of maximal rank, 
the emergence of the striped phase can also be conjectured from the eigenvalue distribution 
of the matrix, which represents the scalar field in NC geometry [18,21]. 

In the case of 4d U(l) NC gauge theories, the IR singularity occurs in the one-loop 
calculation of the vacuum polarization tensor [22]. This changes the dispersion relation of 
photons at low momenta [23]. In fact the IR singularity in the dispersion relation occurs 
with an overall negative sign [24-26] unless one introduces sufficiently many fermionic 
degrees of freedom. It turns out that the manifestly gauge invariant effective action, which 
is derived from field theoretical calculations as well as from string theoretical calculations 
[27,28], involves the open Wilson line operators [29], which implies that the IR singularity is 
associated with their condensation. This causes the spontaneous breakdown of translational 
invariance since the open Wilson lines carry specific non-zero momenta dictated by the 
star-gauge invariance [29]. In the case of a finite NC torus, the same phenomenon has 
been understood [30] from its Morita dual commutative SU{N) gauge theory with twisted 

^See refs. [11,12] for related analytic works, including several attempts to actually solve this model. 
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boundary conditions^. 

We should note, however, that these perturbative analyses merely imply that the per- 
turbative vacuum of 4d NC gauge theory has tachyonic instability, but they do not exclude 
the possibility that the system eventually finds a stable vacuum. (An analogous scenario 
known as the "tachyon condensation" has been proposed in the open string field the- 
ory [31]^.) In order to investigate this issue, we definitely need to employ a non-perturbative 
framework. 

In the current work, the lattice formulation [7] is particularly important since it allows 
us to study the model from first principles by Monte Carlo simulation. We study pure 
U(l) gauge theory in four dimensions, where two directions are non-commutative, and we 
find a phase in which open Wilson lines acquire non-zero vacuum expectation values. In 
fact this phase extends towards weak coupling as the system size is increased, and it turns 
out that we always end up in this phase in the continuum and infinite-volume limits with 
a fixed non-commutativity parameter. This is consistent with the prediction from the one- 
loop calculations. We observe, however, that various physical quantities obey clear scaling 
behaviors in the above limit, which provides evidence for a possible continuum theory. We 
also study the dispersion relation. In the weak coupling symmetric phase, it involves a 
negative IR-singular term, which makes the energy vanish at a non-zero momentum as one 
approaches the critical point. In the broken phase, although the momentum components 
in the NC directions are no longer conserved, we can still study the relation between the 
energy and the momentum in the commutative direction. This reveals the existence of the 
Nambu-Goldstone mode associated with the spontaneous symmetry breakdown. 

This paper is organized as follows. In section ^ we define a gauge theory in 4d NC space, 
and we briefly review the results obtained by perturbative calculations. In section ^ we 
introduce the lattice formulation and discuss how to take the continuum limit. In section ^ 
we present the phase diagram of the lattice model. In section ^ we investigate the existence 
of a sensible continuum limit. In section ^ we study how the space in the NC directions 
looks like in the broken phase. In section ^ we study the dispersion relation in each 
phase to gain deeper understanding of the observed phase transition and the continuum 
limit. Section ^ is devoted to a summary and discussions. In the appendix we explain the 
algorithm used for our Monte Carlo simulations. Some preliminary results of this work 
have been presented in proceeding contributions [34]. 

2. Tachyonic instability in perturbation theory 

NC geometry is characterized by the following commutation relation among the space-time 
coordinates 

[x^„Xu] =i&f_iu , (2.1) 

^At high tempertature (on the scale of the non-commutativity) the ring diagrams dominate. Their 
evaluation to all orders implies a critical temperature, above which the magnetic photon mass squared 
turns negative in the NC planes [32]. 

^This might be more than a plain analogy since the tachyonic instability of NC gauge theory has partial 
contribution from the closed string tachyon [28,33]. 
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where Q^y is the non-commutativity tensor. Here we consider the 4d Euchdean space-time 
(n = 1, • ■ ■ ) 4) with the non-commutativity introduced only in the /i = 1, 2 directions, 

i.e., 

912 = -621 = , 

0/^1/ = otherwise . (2-2) 

Since we have two commutative directions, we may regard one of them as the Euchdean 
time. This allows us to alleviate the well-known problems concerning causality and unitar- 
ity [35,36]. Moreover, it also enables us to define the dispersion relation as a useful probe 
to study the system [19]. 

Pure U(l) gauge theory in this NC space can be formulated in the path integral 
formalism by the gauge action 

F^^ = d^Au - duA^ + ig{A^^ Au- Au^ A^) , (2.3) 

where the space-time coordinates x^ are c-numbers as in ordinary field theories, but the 
non-commutativity is now encoded in the star product 



f{x)^g{x) = exp (^e,.^ A) 



(2.4) 



x=y 



The action in eq. (^^) is invariant under a star-gauge transformation 

A^^ A^ + d^K + ig[A^,K\^ , (2.5) 

where the star-commutator is defined by [/, 5]* = f * g — g * f ■ Note that non-linear terms 
appear in the field strength tensor in eq. ( |2.3| ) due to the NC geometry, although we are 
dealing with a gauge group of rank 1. As a result, the theory shares such properties^ as 
a negative beta function [22,25] with non-Abelian (rather than Abelian) gauge theories in 
the commutative space. 

This model has been studied extensively in perturbation theory. In particular the 
result that is most relevant for us is the one-loop effective action for the gauge field, which 
involves the quadratic term [22-26] 

where = Q^uPu ■ 

This term emerges from the non-planar diagrams and represents an effect of UV/IR mixing. 
Since the effective potential, which is minus the effective action according to the adopted 
convention, involves a negative quadratic term, we find that the low momentum modes of 
the gauge field cause a tachyonic instability. 



*It is also conjectured that NC photons form bound states analogous to glueballs in QCD [37]. 
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In fact the quadratic term in ( p. 61) is invariant under the star-gauge transformation 
( |2.5| ) only at the leading order in A^. In order to obtain a fully gauge invariant effective 
action, one needs to include terms at higher orders in Af^ as in the supersymmetric case [38]. 
The term in the gauge-invariant effective action which is most singular at small p is given 
by [27,28]5 

ri-ioop = ^/J^H^(pm-p)^. (2.7) 

The open Wilson line W{p), which appears here, is a manifestly star-gauge invariant op- 
erator defined as [29] 

W{p) = I d^x e^f-^ V exp, (^i g J^^'^ A^{0 ^^Cm) > (2-8) 

where V exp represents the path-ordered exponential, and the path for the line integral 
over ^ is taken to be a straight line connecting x and x + p. Expanding ( |2.7| ) in terms of 
A^, one obtains ( |2.6| ) at the leading order, up to an irrelevant constant term. The result 
( |2.7| ) implies that the instability is associated with the condensation of open Wilson lines 
W{p) with small p. Since the open Wilson line W{p) carries non-zero momentum, its 
condensation causes the spontaneous breakdown of translational symmetry. 

From the perturbative analysis alone, one cannot tell whether the theory possesses a 
stable non-perturbative vacuum after the condensation of open Wilson lines. To answer 
this question and to study the nature of the theory at the stable vacuum, if it exists, we 
definitely need a fully non-perturbative approach. 

3. The lattice model and its continuum limit 
3.1 Lattice regularization of NC gauge theory 

The lattice regularized version of the theory (|2.3D can be defined by an analog of Wilson's 



plaquette action [7] 

S = -P^^Ufj,{x) -kU^{x + afi) *Uf,{x + au)* -kU,y{x)* + c.c. , (3.1) 

X fl<U 

where the symbol ft represents a unit vector in the /x-direction and we have introduced the 
lattice spacing a. The link variables U^{x) (^ = 1, • • • ,4) are complex fields on the lattice 
satisfying the star-unitarity condition 

u^{x) * u^{xy = u^{xy * u^{x) = i . (3.2) 



The star product on the lattice can be obtained by rewriting (2^) in terms of Fourier 
modes and restricting the momenta to the Brillouin zone. The action (3.1) is invariant 
under a star-gauge transformation 

U^{x) ^ g{x) -k U^{x) -k g{x + aft)* , (3.3) 



^Refs. [39] pointed out that the appearance of open Wilson hnes in the effective action is a generic 
feature of field theories on NC geometry. 
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where also g{x) obeys the star-unitarity condition 

g{x)*g{xy =g{xr*g{x) = l . (3.4) 

As in the commutative space, one obtains the continuum action ( |2.3| ) from ( |3.1| ) in the 
a — > limit with the identification d = ttv and 

(r-x+aji \ 
ig J A^iOd^J . (3.5) 

In order to study the lattice NC theory (3.1) by Monte Carlo simulations, it is cru- 
cial to reformulate it in terms of matrices [7]. In the present setup ( |2.2| ), with two NC 
directions and two commutative ones, the transcription applies only to the NC directions, 
and the commutative directions remain untouched. Let us decompose the four-dimensional 
coordinate as x = {y,z), where 

y = {xi,X2) and z = {x3,X4) (3.6) 

represent two-dimensional coordinates in the NC and in the commutative plane, respec- 
tively. We use a one-to-one map between a field ip{x) on the four-dimensional N x N x Lx L 
lattice and a N x N matrix field ip{z) on a two-dimensional L x L lattice. This map yields 
the following correspondence 

ipi{y,z)'kip2{y,z) ipi{z)ip2{z) , (3.7) 

ip{y + a[i,z) T^tp{z)Tl^ , (3.8) 

l^^^iy,z) ^ ^tr^iz). (3.9) 

y 

The SU(A^) matrices (/i = 1, 2), which represent a shift in a NC direction, satisfy the 't 
Hooft-Weyl algebra 

riT2 = Z12T2T1 , (3.10) 
2i2 = = exp ^vri ^ ^ , (3.11) 

with the matrix size being an odd integer. For this particular construction [10,19,40], 
which we are going to use throughout this paper, it turns out that the non-commutativity 
parameter in (^]^) is given by 

e = -Na'^ . (3.12) 
vr 

Note that the extent in the NC directions Na goes to infinity in the continuum limit a — > 
at fixed 9. 

Using this map, the link variables U^{x) are mapped to a x matrix field U^{z) on 
the two-dimensional L x L lattice, and the star-unitarity condition (^]^) simply requires 
U^{z) to be unitary. The action (3T) can be rewritten in terms of the unitary matrix field 
Ufi{z) in a straightforward manner. By performing a field redefinition 

U^{z) for // = 3,4 , 
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where V^{z) G U(A^), we arrive at 



•S" = 5'nc + •5*00111 ~l~ 'S'mixed > 
Snc = -NpZuY^tv (Vi{z)V2{z)Vi{z)W2{zy)+ c.c. , 



iV/3 ^ tr ( Vsiz) V^iz + a3) V^iz + a4)t y4(z)t) + c.c 



z 

2 4 



'mixed 



-NP^^^tT [V^{z)V,{z)V^{z + ai)yV,{z)^j + c.c. . (3.14) 

z ^=1 j/=3 



These terms represent a plaquette in the NC plane, in the commutative plane and a sum 
over the plaquettes in the four mixed planes, respectively. Since the action ( 3.14| ) is invari- 
ant under the gauge transformation 

[aizjyf^iz) 9{z + afj.y for /i = 3,4, 

where g{z) G U(A^), it defines a 2d lattice gauge theory, in which V^{z) (fi = 3,4) cor- 
respond to link variables, and V^iz) {fi = 1,2) correspond to scalar fields in the adjoint 
representation with the specific self-coupling given by 5'nc- Taking into account the field 
redefinition ( |3.13| ), one easily finds that the U(A^) gauge symmetry (3.15) of the 2d theory 



corresponds precisely to the star-gauge invariance ( |3.3| ) of the 4d NC theory that we started 
with. 

3.2 Eguchi-Kawai equivalence in the planar limit 



In fact the lattice model ( 3.14 ) can be obtained by the twisted dimensional reduction [8] 
from 4d pure U(A'^) lattice gauge theory^. Historically such a model appeared in the context 
of the Eguchi-Kawai reduction [42] , which in the present case provides an explicit relation 
between the 2d U(A'^) theory ( 3.14| ) and the 4d U(A^) theory (both in the commutative 



space) in the large- limit at fixed /3. (This limit is usually referred to as the planar limit 
since only planar diagrams survive [43].) More specifically, the vacuum expectation value 
of a Wilson loop defined in each theory coincides in the above limit, if the U(l)^ symmetry 

V^{z) ^ e^'^V^iz) for ^ = 1, 2 (3.16) 



of the 2d theory (3.14) is not spontaneously broken. We will see that this condition holds 
in the weak coupling and the strong coupling regions, but it is violated in the intermediate 
region. In fact this symmetry ( p. 16 ) includes the (discrete) translational symmetry of the 



4d NC lattice gauge theory in the NC directions as a subgroup up to a star-gauge trans- 
formation [44]. Therefore, the spontaneous breaking of the U(l)^ symmetry corresponds 
precisely to the IR instability of the perturbative vacuum discussed in section |2[ 



®In the present case, this amounts to considering pure U(A^) gauge theory onalxlxLxL lattice with 
the twisted boundary condition [41] for the 1,2 directions, and the periodic boundary condition for the 3,4 
directions. The phase factor Z12 in the first term of (p.l^) represents the twist. 
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Let us comment on some known results for a totally reduced (one-site) model. The 
condition for the Eguchi-Kawai equivalence in this case is that the U(l)^ symmetry is not 
spontaneously broken. The twist was introduced in ref. [8] into the original (untwisted) 
Eguchi-Kawai model [42] in order to avoid the spontaneous U(l)^ symmetry breaking in 
the weak coupling region^. Indeed the condition is satisfied in the strong coupling and 
weak coupling expansions, and early Monte Carlo studies suggested that it holds, too, at 
intermediate couplings [8]. However, Ishikawa and Okawa [47] have recently performed 
Monte Carlo simulations with much larger N, and observed a signal of spontaneous sym- 
metry breaking at intermediate couplings. This behavior persisted also for choices of the 
twists other than the minimal one adopted in ref. [8]. 

The Eguchi-Kawai equivalence, if the condition is met, holds in the large- limit with 
a finite lattice spacing a. The continuum limit a — > as a commutative 4d U(oo) gauge 
theory can be taken in the next step by sending /3 — > oo. In our model, however, the critical 
P, below which the symmetry breaking occurs, is observed to increase as A^^, which implies 
that the Eguchi-Kawai equivalence does not persist in the continuum limit. It remains 
to be seen whether this is the case even if one takes other options in defining twisted 
reduced models (the choice of the twist, partial or total reduction, etc.). An alternative to 
ensure the validity of the Eguchi-Kawai equivalence in the continuum limit is the quenched 
reduced model [48]. Its continuum version [49] has been applied recently [50] to a non- 
lattice regularization of Matrix Theory [51]. See also refs. [45,52] for recent developments 
in the large- reduction in the continuum. 

3.3 Double scaling limit 

As one can see from ( [3.12 ), the planar limit corresponds to = oo in the context of NC field 



theory. In order to take the continuum limit as a NC gauge theory with finite 9, one has to 



take the — > oo and /3 — > oo limits simultaneously, while satisfying ( 3.12 ), which therefore 
implies a ^ 0. Following the usual terminology in matrix models, we call it the double 
scaling limit. Unlike in the planar limit, non-planar diagrams survive and they cause the 
intriguing UV/IR mixing effects. The inverse coupling constant /? has to be tuned in such 
a way that physical quantities converge. Whether this is possible or not is precisely the 
issue of (non-perturbative) renormalizability, which we address in the following. 

In the 2d NC gauge theory [10], since the U(l)^ symmetry is not spontaneously broken 
throughout the whole coupling region, the Eguchi-Kawai equivalence holds in the planar 
limit for all /?. The tuning of /? after taking the planar limit can therefore be deduced from 
the exact solution [53] of the 2d U(cxd) lattice gauge theory. This tuning of P with respect 
to the lattice spacing a should be used also in the double scaling limit in order to make 
the correlation functions scale in the UV regime [54]. Whether the scaling extends to the 
IR regime or not is a non-trivial issue, which was answered affirmatively in ref. [10]. In the 
present 4d case, since the U(l)'^ symmetry is spontaneously broken in the coupling region 



^It has been shown recently [46] that in the original (untwisted) Eguchi-Kawai model and in a 
more general model, the U(l)* symmetry breaking does not occur at once, but rather in a sequence 
U(l)*— >U(1)^— >U(1)^— >U(1) ^none. Such a partial breaking of U(l)* symmetry was observed earlier 
in large A'^ gauge theories in a finite box [45]. 
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relevant to the continuum limit, the tuning of f3 cannot be deduced from the results in 4d 
U(oo) lattice gauge theory. We therefore have to fine-tune /? in such a way that the double 
scaling is optimized. 

4. Phase structure 



To begin with, we investigate the phase structure of the lattice model (3.14). Throughout 
this paper, we take L, the number of sites in the commutative directions, to be L = it 1 
so that it becomes a multiple of four in order to use four processors on a parallel machine 
efficiently. The small anisotropy of the lattice should be negligible at large A^. 

As a standard quantity we plot the normalized action (or the average plaquette) 

1 



S 



-AS) 



(4.1) 



12A^2^ 

against P for A^ = 25 in figure |l|. We have performed a "thermal cycle" by carrying 
out simulations at varying (3 with the initial configuration taken from a configuration 
thermalized with a slightly larger /smaller /3. There is a gap at /3 ~ 0.35 (left), and we 
observe a slight tendency of a possible hysteresis behavior at 1-1 < /? < 1.6 (right) although 
the difference between the two branches, even it exists, is too tiny to be confirmed definitely. 
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Figure 1: (Left) The expectation value of the normalized action (4.1) is plotted against P for 
A^ = 25. The two kinds of symbols are used to distinguish the results obtained for decreasing f3 
and increasing (3 in the "thermal cycle" . The solid lines represent the results of the strong and the 
weak coupling expansions. The discontinuity is observed at /3 ~ 0.35. (Right) The zoom-up of the 
same plot around 1.1 < /3 < 1.6. 



As an order parameter for the spontaneous breaking of the U(l)2 symmetry ( p. 16 ), we 
define a gauge invariant operator 

^ ^trfy^(z)") for /x = 1,2, (4.2) 



which transforms non-trivially under ( |3.16| ). This operator corresponds to the open Wilson 
line (]2.8| ) carrying a momentum with the absolute value [7, 10] 



^ for even n , 
^ for odd n . 



(4.3) 
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Since the operator P^{n) with odd n carries a momentum of the cutoff order, it does not 
couple to excitations that survive in the continuum Umit^. Therefore we wih focus mainly 
on the even n case in what follows. 
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Figure 2: The order parameter is plotted against (3 for n = 2 (left) and n — 4 (right). 

The system size is = 15 (circles), N — 25 (triangles) and = 35 (squares). The closed (open) 
symbols represent results obtained with increasing (decreasing) /3, which show a clear hysteresis 
behavior. 



In figure § we plot (|Pi(n)|) against 
/3 for n = 2 (left) and n = 4 (right). 
Since it turns out that |Pi(n)| and |P2(?t-)| 
are not too different for each configura- 
tion^, we actually take an average over 
the two NC directions to increase the 
statistics. We observe that there is a 
phase, in which the order parameter be- 
comes non-zero. On the other hand, the 
quantity (|Pi(n)|) for odd n takes tiny 
values throughout the whole region of 
/?. This implies that the U(l)^ sym- 
metry is broken down to (Z2)^ in this 
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Figure 3: The upper (squares) and lower (circles) 
critical points between the symmetric phase and the 
broken phase, obtained from figure 0. The lines re- 



phase, which we refer to as the "broken ^^^^^^^ a fit to /3 = ciiV^ + cs, where ci = 0.00226(1), 



C2 — 0.145(4) for the upper critical point, and ci = 
0.00141(1), C2 = 0.250(5) for the lower critical point. 



phase". The fact that the U(l) sym- 
metry is unbroken, both in the small (5 
regime and the large /3 regime, can be 

*This 7i2 grading of the open Wilson line operators should be regarded as an artifact of the lattice 
regularization [7]. 

^In principle, it could also be that the U(l)^ symmetry is broken down only to U(l). In that case |Pi(n)| 
and |P2(?i)| would be very different for each configuration, and taking the average over the two directions 
would not be legitimate. We observed that this is not the case. 

^"A similar phase structure is obtained in the totally reduced model with the minimal twist [47 
that case, however. 



In 



Pi(n)l) for odd n also acquires a non-zero value. This is not so surprising since the 
argument based on the momentum spectrum (4.3) does not apply to the model with the minimal twist. 
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understood from the strong coupling and the weak couphng expansions [8] as in the totally 
reduced model^''. 

The transition between the strong coupling phase and the broken phase occurs at 
(3 ~ 0.35 (for all N), which coincides with the position of the gap in figure |l| (left). Actually 
this value agrees with the critical point of the bulk transition known in the 4d SU(A^) lattice 
gauge theory with Wilson's plaquette action at large [55], which is also reproduced by 
large- A^ reduced models [8,56]. Since the U(l)^ symmetry is unbroken in our model in 
the strong coupling phase, something must occur at the critical /3 of the bulk transition 
from the viewpoint of the Eguchi-Kawai equivalence. What actually happens is that the 
U(l)^ symmetry is spontaneously broken down to (22)^, and the equivalence ceases to hold 
precisely at that point. 

The transition between the broken phase and the weak coupling phase is more impor- 
tant since it is relevant for the continuum limit. The hysteresis behavior, which indicates 
a first order phase transition, can now be seen clearly, unlike in figure ||. In fact the lower 
critical point of this phase transition, denoted as /3c in the following, can be estimated by 
perturbation theory, and we obtain [57] j3c ~ A^ at large A, which is consistent with our 
results shown in figure ^. This implies, in particular, that one always ends up in the broken 
phase if one takes the A^ — > 00 limit at fixed (3 (planar limit), which corresponds to = cxd 
in the context of NC gauge theory. 

Let us next consider closed Wilson loops, which play an important role in commutative 
gauge theories as a criterion for confinement. In the present case, since we introduce non- 
commutativity only in two directions, there are three kinds of square-shaped Wilson loops 
depending on their orientations. Using the parallel transporter 

V^(z,n) = V^{z)Vy{z + ai))---V^{z + {n-l)ai)) (4.4) 

in the commutative directions, the closed, square-shaped Wilson loops can be defined as 

z 

W^u{n) = ^^tr(y^(z)"V.(z,n)y^(z + naz))t"V.(z,n)t) , 

z 

Wu{n) = ^^tr (v3(z,n)V4(z + na3,n)V3(z + na4,n)tV4(z,n)t) , (4.5) 
2 

where /i = 1, 2 and = 3,4. We may define W^y{n) for /i > z/ in a similar manner, but it 
suffices to consider ( ^^ ) since Wf^^^n) = Wu^{n)*. Mapping the expression ( [4.5| ) back to 
the 4d NC space, we find that the operator W^^in) represents a closed n x n Wilson loop, 
whose starting point is integrated over the 4d space. Therefore it does not carry non-zero 
momentum unlike the open Wilson line P^{n). This corresponds to the fact that W^^in) 
is invariant under the U(l)^ transformation ( |3.16| ), whereas P^{n) is not. 

Let us consider the expectation value (VF^,y(n)). First we reduce the number of obser- 
vables using discrete symmetries. Due to the symmetry under exchanging the third and 
the fourth directions, we have (W34) = (W43) G M and (W^s) = (W^4), where = 1,2. 
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Figure 4: (Left column) The expectation value of the Wilson loop is plotted against its size n 
for various P at N — 35. The results for f3 = 2.0 are obtained in the symmetric phase. (Right 
column) The expectation value of the Wilson loop is plotted against its size n for N = 15, 25, 35, 45 
a.t /3 — 1.5. The results for N — 15, 25 are obtained in the symmetric phase. 



Based on the symmetry under Vi{z) i-^ V2{z), V2{z) ^ Vi{z)\ we obtain {Wiu) = {W2u) £ 
M, where v = 3,4. It then follows that {W^u) with fi = 1,2 and v = 3,4 are all real 
and equal. We will use (Wi^) as a representative, but in actual simulation we measure 
I ^^^3(Re VF^i,) to increase the statistics. Note that {W12), which represents a 

closed Wilson loop in the NC directions, is complex in general unlike in commutative 
gauge theories. This occurs because the non-commutativity tensor Q^u breaks parity. 

In figure ^ (left column) the expectation value of the Wilson loop is plotted against its 
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linear size n for various (3 values at = 35. First let us focus on the results at (5 = 2.0, 
which are obtained in the symmetric phase. Here the three types of Wilson loop are almost 
identical and follow the perimeter law over the whole range of n shown in the figure. This 
implies that the system is qualitatively similar to 4d compact U(l) lattice gauge theory in 
the commutative space, which is non-confining at weak coupling. Note, however, that the 



effect of the NC geometry is seen in the dispersion relation for the same (3 (figure 11). As 
we decrease /3, the system enters the broken phase, and the three kinds of Wilson loops 
start to drift apart. 

The right column of figure § shows the results for various N at /3 = 1.5. The data 
for = 15 and N = 25, which are obtained in the symmetric phase, lie almost on top of 
each other. This can be interpreted as the scaling behavior corresponding to the would-be 
planar limit in the symmetric phase, although one actually enters the broken phase as A'^ is 
increased further. In fact the effect of increasing A'^ is similar to that of decreasing (5, which 
suggests the possibility to make the results scale by increasing (3 and A^ simultaneously. 
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Figure 5: The expectation value of the Wilson loop 
in the commutative plane. By tuning (3 depending 
on iV, we can make the results scale. 



5. Existence of a continuum limit 



In this section we investigate whether 
it is possible to tune (3 function of 
A^ in such a way that physical quantities 
scale. For the rest of this paper except 
for figure |^, we set a = 1 for A^ = 45 
as a convention, and the lattice spacing 
a for other A^ is determined^ ^ through 
( pT2D with 6^ = f =i 14.3. 

As a practical strategy to fine-tune 
/3, we focus on the expectation value of 
the square Wilson loop ^¥34(71) in the 
commutative directions since it has a 
smooth dependence on its size n. Ta- 
ble |l| shows the optimal values of j3 for 
each A^, and figure ^ shows the corresponding plot. The horizontal axis represents the 
physical size (no) of the loop. We observe a clear scaling behavior, and the scaling region 
extends as A^ increases. In fact the optimal (3 increases with A^ much slower than the lower 
critical point [3c between the broken phase and the weak coupling phase, which grows as 
A^^ (see figure ^). This implies that we remain in the broken phase in the double scaling 
limit. 

Let us see whether other quantities scale as well using the same sets of parameters 
as those given in table |l[ In figure |^ we plot the expectation value of the Wilson loop in 
the NC plane (left) and in the mixed planes (right). We do observe a compelling scaling 
behavior. 

The Wilson loop in the commutative plane follows the perimeter law at large size, which 
suggests that the theory is non-confi ning in the commutative directions. The Wilson loop 
^^If we had discarded the constraint ( 3.12 ) and determined a for each A'^ to optimize the scaling, we 
would have concluded that 6 had to be fixed in order to obtain the scaling behavior. This follows from the 
fact that the scaling functions obtained for fixed 9 have non-trivial ^-dependence. That 6 does not receive 
renormalization should therefore be considered as an observation rather than an input of our study. 



in the NC plane is complex in general, and its real part oscillates around zero. In figure 
1^ we plot the absolute value and the phase of the Wilson loop against the physical area 
A = {na)'^. The absolute value decreases monotonously obeying roughly the area law. This 
is not so surprising since non-Abelian nature comes in through the star product in 
although we are dealing with a U(l) theory. The phase <I> grows linearly as <I> = ^A, which 
is reminiscent of the Aharonov-Bohm effect with the magnetic field B = 6^^ [4]. The same 
behavior has been observed in the 2d case [10]. 

The behavior of the Wilson loop in the mixed directions 
is somehow between the other two kinds. It decreases fol- 
lowing roughly the perimeter law for large size, but a slight 
oscillating behavior seems to be superimposed (although these 
Wilson loops are real). This suggests that the shape of the 
potential between a quark and an anti-quark separated in a 
NC direction is oscillating, but we do not have a clear physical 
interpretation of such a behavior. Table 1: The sets of param- 

Let us turn to the open Wilson line, which was studied in eters used for the double seal- 
section ^ to investigate the spontaneous breakdown of U(l)^ i^S limit, 
symmetry. In figure |8| (left) we plot the expectation value 

(|Pi(n)|) for even n against the momentum defined by ( [4.3[ ). We observe a tendency that 
the results lie on a single curve except for the n = 2 data (corresponding to the point at 
the smallest p for each A^). A similar anomalous behavior is seen also in figure ^ for n = 3 
data. We consider these behaviors to be finite A^ artifacts since they tend to disappear 
with increasing A^. 
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Figure 6: The expeetation value of the Wilson loop in the NC plane (left) and in the mixed planes 
(right). The values of /3 are given in table which was chosen in such a way that the expectation 
value of the Wilson loop in the commutative plane scales. 



6. Extent of the dynamical NC space 

Since the translational invariance in the NC directions is spontaneously broken in the 
broken phase, it is natural to ask how the space actually looks like in those directions. 
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Figure 7: The absolute value (left) and the phase (right) of the Wilson loop in the NC plane 
is plotted against the physical area A = (na)^. The absolute value follows roughly an area law. 
Beyond small areas, the phase $ agrees with the Aharonov-Bohm like behavior ^ — A/9 represented 
by the solid straight line. 



As a related quantity, let us consider the eigenvalues of the unitary matrix V^{z) 
(n = 1,2), which we denote as e*'^™(^) (j = 1, • • • , A^), where — vr < '&^j{z) < vr. The open 



Wilson line (4.2) can be written in terms of the eigenvalues as 



1 ^ 

^/^(") = ]^EE«^"'"^^^- (6-1) 

z j=l 



Since the U(l) transformation (|3.16D rotates all the eigenvalues by a constant angle a^, 
the eigenvalue distribution should be uniform if the U(l)^ symmetry is not spontaneously 
broken. This is how the open Wilson line serves as an order parameter. In the broken 
phase, we have seen that the open Wilson line Pfi{n) acquires a non-zero expectation value 
only for even n. This implies that the eigenvalues are clustered in two bunches in a Z2- 
symmetric way. Since the translation in the NC direction is represented by the U(l)'^ 
transformation, we may consider that the eigenvalue distribution represents the "shape" 
of the NC space. 

Let us define the eigenvalue distribution by 

z ^=1 j=l ^ ^ ^ ' 

where we have taken an average over the two NC directions as we did for (|Pi(n)|). The 
coefficient of 'd^j{z) is motivated from the corresponding relation p = -^n {n : even) for 
the momentum conjugate to the coordinate x. Before taking an ensemble average, we 
rotate each configuration according to (|3.16| ) in such a way that X^^tr V^(z)^ becomes real 
positive^^. Since the eigenvalue distribution p{x) is invariant under the shift x ^ x + Na 
modulo 2A^a, we may restrict ourselves to the fundamental domain \x\ < Na/2. Figure 
|8| (right) demonstrates a clear scaling behavior of p{x). The eigenvalue distribution can 
be interpreted as the density distribution of D-branes [27]. If we take the view point of 

^■^We average over the Z2-ambiguity in fixing the angle. 
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Figure 8: (Left) The expectation value (|-P^(«)|) for even n is plotted against the momentum 
p = -^n for the parameters listed in table |l|. (Right) The eigenvalue distribution p{x) is plotted 
for the parameters listed in table 0. 



string theory, in which the space-time is represented by the eigenvalue distribution of the 
matrices [61,62], our results imply that the "dynamical space" in the NC directions has 
shrunk, but it has a finite extent in the double scaling limit. 

Since we have seen that the extent of the "dynamical space" in the NC direction is 
finite for a fixed 9, it is natural to ask how the extent depends on 9. Let us note that 

1 2 i.Na/2 

-^{\P,{2)\) = 2 dxe^^pix) 

2 ^ J-Na/2 

where taking the absolute value on the left-hand side corresponds to rotating 'd^j{z) in 
the definition ( |6.2D of p{x) before taking the ensemble average as explained above. In the 
last line, we have assumed that p{x) is peaked around x = in the fundamental domain 
|x| < Na/2. 

As a definition for the extent of the dynamical space in the NC direction, we therefore 

use 




It should be mentioned that the space-time, on which the NC gauge theory is defined, has 
the extent Na, which diverges in the double scaling limit. However, since the translational 
invariance in the NC directions is spontaneously broken, the observer on the NC space-time 
does not recognize that the space-time extends to infinity. The quantity Ax measures the 
extent of the space-time in the NC directions, which is recognized by the observer to be 
qualitatively uniform. In that sense. Ax is analogous to the width of the stripes [19] in the 
3d NC Xcj)^ theory, in which the translational invariance in the NC directions is broken due 
to the space-dependent vacuum expectation value of the scalar field. 
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Instead of repeating the whole pro- 
cedure of taking the double scaling limit 
at each 0, here we assume that the dou- 
ble scaling is obtained for the same sets^^ 
of (a, (3) listed in table ||. This assump- 
tion is justified at large N given that the 
ultraviolet properties of NC theories are 
independent of 6. For various values of 
(3 in the range 0.92 < /3 < 2.00, we de- 
termine a by interpolating the relation 
between a and f3 presented in table ||. 
The value of 9 is then determined from 




Figure 9: The extent of the dynamical space in the 



NC directions defined by (6.4) from the view point 
of string theory, is plotted against ^/O for N — 35. 



( 3.12 ) using o and N . In figure ^ we plot 
Ax against ^/9 for = 35. At large 9 
the extent A.x increases linearly with as expected on dimensional account. 



7. Dispersion relation 

In this section we investigate the dispersion relation in the symmetric phase and the bro- 
ken phase separately. The analysis in the symmetric phase reveals the IR singularity (^]^) , 
which is responsible for the phase transition, whereas the analysis in the broken phase 
enables us to identify the Nambu-Goldstone mode associated with the spontaneous break- 
down of the U(l)^ symmetry. Thanks to the existence of the commutative directions in the 
present setup, we may regard one of the coordinates (say, XtC) as "time". From the expo- 
nential decay of the two-point correlation function of open Wilson line operators separated 
in the "time" direction, we can extract the energy of a state that couples to the operator. 
Similar studies have been done also in the case of NC scalar field theory [19]. 

7.1 Results in the symmetric phase — IR singularity 

Let us define the open Wilson line operator at a fixed time 3:4 as 



P^(x4,n) = ^^tr 

^3 



V^.{xz,XiT] for ^=1,2 



(7.1) 



which has a zero momentum component in the X3 direction^^ and a non-zero momen- 
tum component (|4.3| ) in a NC direction depending on n. Then we define the two-point 
correlation function of the open Wilson lines 

, 2 

Cn{r) ^-Y. E(^m(^4, nf . P,ix, + r, n)) (7.2) 
[1=1 X4 



^^In fact one can always multiply the lattice spacing a by a ^-dependent factor without affecting the 
scaling property. This ambiguity corresponds to the arbitrary choice of the A-parameter at each 6. The 
assumption we adopted here corresponds to taking the A-parameter independent of 9. 

^*One can introduce a non-zero momentum component ps in the 2:3 direction by inserting a phase factor 
g»P3a;3 jjj ^}jg summation over 3:3 in eq. ([7.l|). The dispersion relation ( [7.3| ) will then have a term (^3)^ on 
the right-hand side. 



-17- 




2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 



Figure 10: Two-point correlation function of the open Wilson lines for even n (left) and odd n 
(right) is plotted against the separation r for = 35 and /3 = 2.00. The vertical axis is normalized 
in such a way that the function starts off from unity at the origin. 



with a separation r in the temporal direction. In actual measurement we also consider 
the case with the roles of X3 and X4 exchanged, and take an average over the two cases to 
increase the statistics. 



In figure 10 we plot the two-point 
correlation function for even n (left) and 
odd n (right). For even n we observe 
clear exponential behaviors e"'*''^, from 
which we extract the energy E = X/a 
at each momentum p = -^n. For odd n 
the decay is very rapid, which suggests 
that the energy (as well as the momen- 
tum) is on the cutoff scale. 
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In figure |llj we show the dispersion 
relation obtained from the two-point cor- 
relation function (^^) for even n. The 
set of parameters {N, a, (3) is chosen as 
in the broken phase. Namely, a is de- 
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Figure 11: The dispersion relation in the symmet- 
ric phase. The energy E obtained from the two- 
point correlation function ( [7.2|) is plotted against 
the momentum p for {N,a,(3) = (35,1.13,2.00), 
(25,1.32,1.69), (15,1.73,1.29). 
14.3, and we fine-tune (3 at each in 

such a way that the scaling behavior of the Wilson loops in the commutative plane is opti- 
mized. It turns out that the data points {E,p) for different N lie to a good approximation 
on a single curve 



P 



{Opf 



(7.3) 



with c ~ 0.1285. This form of the dispersion relation is expected from the one-loop cal- 
culation of the vacuum polarization [23-26] . Due to the negative sign of the second term 



in (O), the usual Lorentz invariant (massless) dispersion relation is bent down. The IR 
singularity is regularized on the finite lattice since the smallest non-zero momentum (which 
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corresponds to n = 2) is given by jj^ oc However, if one attempts to increase N fur- 
ther, the energy at the smallest non-zero momentum vanishes at some A^, and one enters 
the broken phase. Therefore we cannot take the double scaling limit in the symmetric 
phase. The scaling observed in the symmetric phase represents an effective theory with 
a finite cutoff. In the case of NC scalar field theory [19], the double scaling limit can be 
taken in the symmetric phase since the IR singularity appears in the dispersion relation 
with the positive sign. 

7.2 Results in the broken phase — Nambu-Goldstone mode 

When we extend the study of the 
dispersion relation to the broken phase, 
we should note that the momentum com- 
ponents in the NC directions are no longer 
conserved, although the momentum com- 
ponent in the commutative direction still 
is. 

Let us therefore define the open Wil- 
son line operator at a fixed time X4 car- 
rying momentum p in the commutative 
direction as 
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Figure 12: The dispersion relation in the broken 
phase. The energy E obtained from the two-point 
correlation function (^^) is plotted against the mo- 
mentum p for the parameters listed in table |^. 



for /i = 1,2. We have chosen the power 

of to be the smallest even number so that the operator couples most effectively to the 
Nambu-Goldstone mode associated with the spontaneous breakdown of the U(l)^ symme- 
try. Then we define the two-point correlation function of the open Wilson lines 



Cp{r) 



^ 2 ^ 



(7.5) 



with a separation r in the temporal direction. As in the symmetric phase, we also consider 
the case with the roles of X3 and X4 exchanged, and take an average over the two cases to 
increase the statistics. 

We measure the two-point correlation function and extract the energy at each mo- 



mentum p. The result is shown in figure 12. It is consistent with the expected massless 
behavior E = p. The discrepancies observed at both ends of the spectrum may be inter- 
preted as finite volume effects and finite lattice-spacing effects, respectively, and they tend 
to disappear as N increases. 

In the 3d NC A(/>^ theory studied in ref. [19], we only have pseudo Nambu-Goldstone 
modes, since the translational symmetry (which is spontaneously broken) is discretized on 
the lattice. The corresponding energy therefore vanishes only in the continuum limit. In 
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the present case of gauge theory, the translational symmetry is enhanced to the continuous 
U(l)^ symmetry even on the lattice. Therefore, the energy corresponding to the Nambu- 
Goldstone mode at zero momentum is exactly zero even before taking the continuum limit. 
This motivated us to introduce a non-zero momentum component p in the commutative 
direction, which makes the energy non-zero. (Note also that for p = 0, we need to subtract 
the disconnected part, which amplifies numerical uncertainties.) The Nambu-Goldstone 
mode can still be identified by studying the dispersion relation. 

8. Summary and discussions 

In this paper we studied four-dimensional gauge theory in NC geometry from first prin- 
ciples based on its lattice formulation. In particular we clarified the fate of the tachyonic 
instability encountered in perturbative calculations. The lattice formulation is suited for 
such a study since the IR singularity responsible for the instability is regularized in a star- 
gauge invariant manner, and we can trace the behavior of the system as the regularization 
is removed. This revealed the existence of a first order phase transition associated with the 
spontaneoTis breakdown of the U(l)^ symmetry, which includes the translational symmetry 
in the NC directions as a discrete subgroup. 

In the weak coupling symmetric phase, we studied the dispersion relation and confirmed 
the presence of an IR singularity. This IR singularity prevents us from taking the continuum 
limit in the symmetric phase. While we are able to observe scaling up to some finite N, we 
cannot increase N further since the energy of the lowest momentum mode vanishes, and 
the vacuum becomes unstable. In the broken phase, however, we provided evidence for a 
sensible continuum limit. The phase is characterized by the condensation of open Wilson 
lines, which represent the "tachyon" in the unstable perturbative vacuum. By studying 
the dispersion relation in the broken phase, we confirmed the appearance of the Nambu- 
Goldstone mode associated with the spontaneous breakdown of the U(l)^ symmetry. 

By measuring the eigenvalue distribution of in the NC directions, we find that the 
dynamical space in the NC directions has shrunk, but it has a finite physical extent in 
the continuum limit. An analogous first order phase transition is found in gauge theories 
on fuzzy manifolds [58], where the fuzzy manifolds collapse at sufficiently large couplings. 
The instability in those cases is due to the uniform condensation of a scalar field on the 
fuzzy manifold^^. The phenomenon that the space-time itself becomes a dynamical object 
is characteristic to gauge theories on NC geometry. This should be closely related to the 
dynamical compactification of extra dimensions'^ in string theory [61, 62] based on the 
IKKT matrix model [63] . Let us recall that the NC geometry appears in string theory as a 
result of introducing a background tensor field [4] . Our conclusion therefore suggests that 
the dynamical compactification in string theory may be associated with the spontaneous 
generation of the background tensor field in six dimensions. This is reminiscent of the 
spontaneous magnetization in 4d non-Abelian gauge theory at finite temperature [64] or 
in 3d QED [65]. 

^^The fuzzy manifolds can be stabilized by adding a sufficiently largo mass to the scalar field [59]. 
'^^As a closely related line of research, see refs. [60], which use fuzzy spheres for compactified dimensions. 



-20- 



On the other hand, if we wish to obtain a phenomenologically viable 4d model, we may 
stay in the symmetric phase by keeping the UV cutoff finite and view the NC gauge theory 
as an effective theory of a more fundamental theory. A 0-deformed dispersion relation for 
the photon such as the one displayed in figure 11 should then have implications [66] on 
observational data from blazars (highly active galactic nuclei) [67], which are assumed to 
emit bursts of photons simultaneously, covering a broad range of energy. Such experimental 
efforts will be intensified in the near future. For instance, the Gamma-ray Large Area 
Space Telescope (GLAST) project [68] is scheduled to be launched in September 2007 and 
to monitor gamma rays from 20 MeV up to 1 TeV. In particular, a relative delay of these 
photons depending on the frequency [69] could hint at a NC geometry. 
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A. The algorithm for the Monte Carlo simulation 



In this section we describe the algorithm used to simulate the model ( 3.14[ ). The main 



part of the simulation is performed by the heat-bath algorithm generalizing ref. [70] , where 
V^{z) is updated by multiplying a SU(A^) matrix. In order to implement the integration 
over the U(A^) — instead of SU(A^) — group manifold appropriately, we also need to in- 
clude a Metropolis procedure for updating V^{z) by multiplying a phase factor These 
procedures, described in the following two subsections, define "one sweep" in our simula- 
tion. Measurements have been performed every 50 sweeps. The number of configurations 
used to obtain an ensemble average is 1980, 1940, 420, 696, 830 for = 25, 35, 45, 55, 65, 
respectively, for the sets of parameters in table ||. 

A.l Heat-bath algorithm for multiplying a SU(A) matrix 

Since two terms, 5nc and 5'mixcd! rnake the action (|3.14D non-linear with respect to V^{z) 
for each fj, and z, we cannot apply the heat-bath method [71] as it stands. Note, however, 
that S'nc is nothing but the action for the twisted Eguchi-Kawai model for each z, which 



''For instance, if we used the heat-bath algorithm alone, det V^(2;) would not change during the simu- 



lation. The Metropolis procedure would be unnecessary if the model (3.14) were an S\J{N) gauge theory 
instead of U(A'^). The difference of U(A'^) and SU(A'^) is irrelevant in the planar limit, but not necessarily 
in the double scaling limit. 
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can be linearized following ref. [70]. Namely we introduce an auxiliary field Qi2iz), which 
is a general complex N x N matrix, and consider the action 

S'^c = NPY1 [tr{Ql2(^)^Ql2(^)} 

-2RetT {Qi2{zy(tVi{z)V2{z)+t*V2{z)Vi{z)'^}] , (A.l) 

where t is a square root of Z12. Completing the squares and integrating out the auxiliary 
field Qi2{z), we retrieve Snc- 



We use a similar trick to linearize 5*1111X6(1 in (3.14). Namely we introduce the auxiliary 
fields Qi3{z), Qii{z), Q23{z) and Q24:{z), each of which is a general complex N x N matrix, 
and consider the action 



z 4 



z ^=1 i/=3 

-2Reti {Q^,{z)^(v^{z)VAz) + VAz)V^{z + ai))^^'j . (A.2) 

Completing the squares and integrating out the auxiliary fields, we retrieve S^nixed- 
Since the new action 

S' = S-^Q + Scorn + 'S'uiixed (^-3) 

is linear with respect to V^{z), we can use the heat-bath algorithm to update V^{z) by 
multiplying a matrix in one of the N{N - l)/2 SU(2) subgroups of the SU(iV) [72]. We 
repeat this procedure for all the independent SU(2) subgroups. The auxiliary fields can be 
easily updated by generating gaussian variables and shifting them appropriately depending 
on the V^{z) field. 

A.2 Metropolis algorithm for multiplying a phase factor 

The updating procedure for rotating the phase of V^{z) is implemented using the Metropolis 
algorithm. It turns out that we can increase the acceptance rate by rotating the phase of 
the auxiliary field Q^u{z) in a covariant manner as 



V^{z) ^ e^°^(")y^(z) f K(z) ^ e^°-^(^)K(^) 

Qi2{z) ^ e^"MWQ,2(z) , 1 Q^u{z) ^ e^°'^(-)Q^.(z) , 



where = 1, 2 and = 3, 4. (No sum is taken in ( A. 4] ) even if an index appears in a term 



twice.) At each step, the transformation parameter a^{z) is non-zero for particular /u and 
z, and it is given by a uniform random number within the range [0, ^], where e = 0.1 
was chosen to keep the acceptance rate reasonably high (e.g., 60%). We accept the trial 
configuration with the probability min( l,e-^^'), ^S[- S'q, where S[ and S'q are the 
action evaluated for the trial configuration and the present configuration, respectively. Note 
that AS' comes solely from S'^j^ed (•S'com) when updating Vfj_{z) with /i = 1, 2 (/x = 3, 4) 



thanks to the phase rotation of the auxiliary field ( A.4 ). We repeat this procedure for all 
choices of /i and z. 
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